Physical Biology

With the development of microfluidics, it has now become possible to assess the susceptibility of bacteria to antibiotics at the single-cell level instead of relying on population measurements. Such studies are particularly relevant when the growth of bacterial population in the presence of antibiotics is heterogeneous. Here, we build a model to describe such a case, and apply it to experimental measurements on a small population of E. Coli exposed to ciprofloxacin, a drug which is well known for triggering a bistable response.

This study introduces a statistical mechanical framework for allosteric communication in proteins based on the spanning-tree ensemble of residue contact networks. By representing protein structures as weighted graphs, we identify each spanning tree as a topological microstate. The canonical partition function is evaluated exactly via the determinant of the reduced weighted Kirchhoff (Laplacian) matrix, allowing for the derivation of global thermodynamic functions (including Helmholtz free energy, internal energy, entropy, and heat capacity) without approximation. Allosteric channels between specific residue pairs are defined as sub-ensembles containing unique simple paths. Using the Burton-Pemantle theorem and the Moore-Penrose pseudoinverse of the graph Laplacian, we compute exact path probabilities and channel-specific thermodynamics. This methodology enables a decomposition of channel heat capacity into energetic and topological components and quantifies residue-level allosteric importance through fractional contributions to the channel partition function. The framework was applied to the G12D mutation in KRAS, comparing wild-type (PDB: 6GOD) and mutant (PDB: 6GOF) proteins. Results show that while the mutation minimally affects mean internal energy and entropy, it reduces global heat capacity by 27.3%. This indicates a topological stiffening where the mutant occupies a significantly narrower landscape of spanning-tree configurations. At the channel level, the mutation maintains distributional stability across six functional routes but triggers a substantial internal redistribution of allosteric importance. Specific residues, such as Q61 and F156, shift occupancy by up to 35.5%. These findings suggest that the G12D mutation does not destroy communication pathways but reorganizes internal information traffic to favor a catalytically impaired state. This approach provides a rigorous, parameter-free metric for understanding how point mutations perturb distal protein signaling.

Cell-cell adhesion is a key regulator of cancer invasion. In this work, we extend a pre-existing individual based cancer invasion model by introducing a stochastic representation of N-cadherin-mediated adhesion, where the lifetime of a cell-cell bond depends on the pulling force acting on the bond. Using experimental data, we derive expressions for the mean and standard deviation of N-cadherin bond lifetimes and fit them to Gamma distributions, enabling their treatment as force-dependent random variables. These distributions are then used to modify the diffusion coefficient of mesenchymal cancer cells. The model predicts reduced random motility with increasing adhesion and incorporates a dynamic transition between catch- and slip-bond behaviour. Along with this model for cell motility, we propose a preliminary physical framework, that can be used to model pattern formation as a result of the new adhesion mechanic.

Understanding how cells migrate through confined environments is crucial for elucidating fundamental biological processes, including cancer invasion, immune surveillance, and tissue morphogenesis. The nucleus, as the largest and stiffest cellular organelle, often limits cellular deformability, making it a key factor in migration through narrow pores or highly constrained spaces. In this work, we introduce a geometric surface partial differential equation (GS-PDE) model in which the cell plasma membrane and nuclear envelope are described as evolving energetic closed surfaces governed by force-balance equations. We replicate the results of a biophysical experiment, in which a microfluidic device is used to impose compressive stresses on cells by driving them through narrow microchannels under a controlled pressure gradient. The model is validated by reproducing cell entry into the microchannels. A parametric sensitivity analysis highlights the dominant influence of specific parameters, whose accurate estimation is essential to faithfully capture the experimental setup. We found that surface tension and confinement geometry emerge as key determinants of translocation efficiency. Although tailored to this specific setup for validation purposes, the framework is sufficiently general to be applied to a broad range of cell mechanics scenarios, providing a robust and flexible tool for investigating the interplay between cell mechanics and confinement. It also offers a solid foundation for future extensions integrating more complex biochemical processes such as active confined migration. Author summaryCells often migrate through very narrow spaces in tissues, a process critical for cancer invasion, immune surveillance, and tissue development. In particular, the stiffness of the nucleus, the largest and most rigid organelle, can limit migration through tight pores. In this study, we present a mathematical model describing the motion of a cell and its nucleus through a microchannel during cell translocation, using a geometric formulation based on surface partial differential equations. The model is general and applicable to a variety of scenarios involving confined cell transport. The model is validated by reproducing key experiments on cell translocation through narrow microchannels. The framework incorporates essential surface features, including mechanical responses, bending rigidity, and surface tension. Sensitivity analysis highlights surface tension and channel geometry as the parameters that most strongly influence translocation. Overall, the model provides new insights into the mechanics of confined cell transport, grants access to cellular quantities that are difficult to measure experimentally, such as cell and nucleus areas, perimeters, and stresses, and establishes a foundation for future extensions incorporating more complex biochemical processes.

Gene drives are genetic elements that can rapidly spread through populations, offering potential solutions for controlling disease vectors and pests. In some scenarios, it is necessary to utilize drives that can be confined to only target populations. The success of these threshold-dependent gene drives, which require a minimum local frequency to establish, depends critically on the spatial strategy used for introduction. Here, we use a reaction-diffusion model to systematically identify optimal release patterns that maximize the per-capita efficiency for four distinct gene drive designs as well as use of Wolbachia bacteria, which spread similarly to frequency-dependent gene drives. We find that the most efficient release strategy is highly dynamic, transitioning from a broad "everywhere" release for short timeframes to a "multiple-ring" pattern for intermediate times, and finally to a focused "center" release for longer timeframes. These timeframes depend on the specific type of drive, with more powerful variants transitioning more quickly to center releases. Our results demonstrate that these optimized, variable release strategies can be substantially more effective than simple uniform releases. This study provides a quantitative framework for designing effective gene drive implementations, highlighting that a carefully planned spatial strategy is essential for maximizing impact, making optimal use of available resources.

A morphospace is an abstract space of theoretically possible biological traits, shapes, or property values. It is interesting to explore which parts of a morphospace life occupies, as compared to those parts which could be occupied, but are not. Comparing random and natural non-coding (nc) RNA secondary structures is an established approach to studying morphospace occupation for RNA structures. Most earlier studies have focused on the minimum free energy (MFE) structure, while relatively few have looked at the Boltzmann distribution, describing the ensemble of energetically suboptimal RNA folds. These suboptimal structures may have important roles and functions, and hence should be examined carefully. Here we compare random and natural ncRNA in terms of their Boltzmann distributions, finding that natural RNA tend to have very similar profiles to random RNA, with the main difference being that natural RNA are slightly more energetically stable, except for very short sequences (20 to 30 nucleotides) which tend to be slightly less stable. We infer that natural ncRNA occupy similar parts of the morphospace that random RNA do, indicating that the biophysics of the genotype-phenotype map largely determines the ensemble properties of ncRNA.

Bacterial communities grow as dynamic populations that respond to their environments. A clinically relevant example is the inactivation of beta-lactam antibiotics by intracellular beta-lactamase in E. faecalis resistant strains. In these populations, resistant bacteria act as antibiotic sinks, detoxifying the environment and allowing sensitive bacteria to survive treatment through a cooperative interaction. In this work, we study strongly coupled planktonic and biofilm populations of mixed sensitive-resistant E. faecalis bacteria under antibiotic stress using fluorescent microscopy. The presence of resistant bacteria in the system benefits both resistant and sensitive cells, leading to mixed planktonic and biofilm populations at super-inhibitory drug concentrations. We show that a beta-lactam antibiotic with or without the addition of a beta-lactam inhibitor can lead to a population inversion effect, characterized by a non-monotonic relation between initial and final fractions of resistant bacteria. The effect is observed in both the planktonic and biofilm populations and is modulated by the total initial cell density. A well-mixed model with competition mediated by resource sharing and cooperation from global degradation of toxins predicts the experimentally observed behavior. These observations suggest underlying population-level mechanisms that are largely independent of biofilm spatial structure.

T cell receptors (TCRs) are critical for immune surveillance and successful adaptive immune response against foreign antigens. TCRs drive this key arm of the immune system through recognition of peptide epitopes presented on MHC complexes. However, they are limited due to their stochastic nature and generation via genetic recombination. In silico design of functional TCRs that target defined peptide epitopes would be of considerable utility but has up until now been unsuccessful. Here, we develop an artificial intelligence (AI)-powered approach using a hybrid physics-based simulation and generative AI that successfully engineers TCRs against defined epitopes presented by MHC-I. We use this approach to design TCRs against two cancer antigens, a HERC1 neoantigen and an immunogenic neoepitope in mutant EGFR. We engineer multiple TCRs against the HERC1 neoantigen which activate T cells in response to exposure to peptide-MHC I and kill cancer cells more effectively than a patient-derived TCR. In addition, we used generative AI to design functional TCRs that target the EGFR T790M neoantigen, engineering greater specificity against the mutant sequence. We present an AI-based approach to TCR design with broad utility for efforts to engineer TCRs and for the development of new cell therapies. One sentence summaryArtificial intelligence-based approach enables the directed engineering of functional TCRs with enhanced features that target cancer neoantigens.

Cells must maintain stable protein levels despite the inherently stochastic nature of gene expression, as excessive fluctuations can disrupt cellular function and impair the reliability of decision-making. Regulatory mechanisms, such as negative feedback, buffer protein fluctuations. Yet, it remains unclear how fluctuations are affected by delays that depend on a molecules specific state. Here, we develop a stochastic model in which proteins are produced in bursts as inactive molecules and pass through a series of intermediate steps before becoming active. The duration of such activation delays depends on the current level of active protein, creating a state-dependent feedback loop. Our model provides explicit analytical expressions relating the delay structure and feedback strength to the variability of active protein levels, quantified using the Fano factor, and shows that state-dependent delays can reduce fluctuations below the baseline expected from simple bursty production. Stochastic simulations confirm these predictions, and incorporating negative feedback in burst production further decreases variability while keeping system behavior predictable. These results reveal how temporal and state-dependent regulation stabilizes protein expression, offering guidance for understanding natural cellular control and designing robust synthetic gene circuits.

Mass-conserved reaction-diffusion systems are used as mathematical models for various phenomena such as cell polarity. Numerical simulations of this system present transient dynamics in which multiple stripe patterns converge to spatially monotonic patterns. Previous studies indicated that the transient dynamics are driven by a mass conservation law and by variations in the amount of substance contained in each pattern, which we refer to as "pattern flux". However, it is challenging to mathematically investigate these pattern dynamics. In this study, we introduce a reaction-diffusion compartment model to investigate the pattern dynamics in view of the conservation law and the pattern flux. This model is defined on multiple intervals (compartments), and diffusive couplings are imposed on each boundary of the compartments. Corresponding to the transient dynamics in the original system, we consider the dynamics around stripe patterns in the compartment model. We derive ordinary differential equations describing the pattern dynamics of the compartment model and analyze the existence and stability of equilibria for the reduced ODE with respect to the boundary parameters. For a specific parameter setting, we obtained results consistent with previous studies. Moreover, we present that the stripe patterns in the compartment model are potentially stabilized by changing the parameter, which is not observed in the original system. We expect that the methodology developed in this paper is extendable to various directions, such as membrane-induced pattern control.

The error catastrophe refers to the proliferation of non-functional molecules in conditions where molecular replication has low accuracy, which is likely to correspond to conditions present at the Origin of Life. This error catastrophe can be avoided thanks to transient compartmentalization, provided that the compartments are sufficiently tight to prevent molecular leakage. Typically, transient compartmentalization models assume that the content of the compartments is completely pooled periodically, resulting in the complete loss of the compositional memory of the compartment. Furthermore, previous models that include the possibility of ecological interactions between molecular parasites and replicators within compartments generally do not study the effect of transient compartmentalization on their coevolution. To address both issues, we develop a framework that accounts for the coevolution of molecular replicators and parasites, along with specific compartmentalization dynamics that are transient yet partial, allowing compositional memory to accumulate from one round of compartmentalization to the next. We benchmark our model with a serial dilution experiment that displays complex oscillatory dynamics among four well-characterized RNA replicators. We also perform experiments to quantify the level of mixing in compartments when stronger stirring tends to homogenize their composition. We then model stirring-induced mixing and show how stirring alters the dynamics of compartmentalized replicators. We conclude that compositional memory arising from transient compartmentalization plays a major role in the dynamics of early molecular systems.

Methods to visualize and quantify the molecular responses of cells to local forces exerted at adhesions are crucial to elucidate how physical forces control cellular behavior. Of the many proteins involved in focal adhesions, vinculin plays a key role in mediating force-sensitive processes. Here, we combined optical tweezers and Forster resonance energy transfer (FRET) microscopy to measure the intensity and FRET efficiency of the vinculin tension sensor, VinTS, in response to a force. Fibroblasts expressing VinTS formed adhesions on fibronectin-coated, 3m-diameter, polystyrene beads. As the beads were displaced by the cell, we applied an optical trap to counteract this movement and increase the traction force required by the cell to maintain the bead displacement. The optical trap stiffness varied from zero (no laser) up to 0.26 pN/nm. In this range, the median bead displacement after 5 min was ~200nm in all trapping conditions inducing counteracting forces in the 10-100pN range. To maintain this displacement, vinculin recruitment increased (up to 35% in relative intensity at high stiffness) while tension increased but more moderately (1-2% decrease in absolute FRET efficiency). For higher trap stiffness, the main response was an increase in vinculin recruitment, while the tension did not increase significantly. The increase in vinculin intensity was correlated with the decrease in FRET efficiency at 0.26 pN/nm but not at lower stiffness. Thus, the presence of the high stiffness optical trap over 5 min appears to induce a positive correlation between vinculin recruitment and vinculin tension. In a few instances, vinculin puncta migrated a few microns away from the bead exceeding the bead movement speed while experiencing an increase in both vinculin intensity and tension. Taken together, the results suggest that combining an optical trap with vinculin tension measurements uncovers novel vinculin dynamics in the presence of a force.

Biological tissues grow at rates that depend on the geometry of the supporting tissue substrate. In this study, we present a novel discrete mathematical model for simulating biological tissue growth in a range of geometries. The discrete model is deterministic and tracks the evolution of the tissue interface by representing it as a chain of individual cells that interact mechanically and simultaneously generate new tissue material. To describe the collective behaviour of cells, we derive a continuum limit description of the discrete model leading to a reaction-diffusion partial differential equation governing the evolution of cell density along the evolving interface. In the continuum limit, the mechanical properties of discrete cells are directly linked to their collective diffusivity, and spatial constraints introduce curvature dependence that is not explicitly incorporated in the discrete model. Numerical simulations of both the discrete and continuum models reproduce the smoothing behaviour observed experimentally with minimal discrepancies between the models. The discrete model offers further individual-level details, including cell trajectory data, for any restoring force law and initial geometry. Where applicable, we discuss how the discrete model and its continuum description can be used to interpret existing experimental observations.

1Correlations between cellular variables, such as gene-expression levels, provide insights into regulatory mechanisms. We focus here on correlations between mRNA and protein levels and re-examine previously derived analytical predictions. We test this prediction on single-cell E. coli data and see substantial disagreement. We hypothesize that this discrepancy arises from the assumption of constant cell volume and develop a theoretical framework for mRNA-protein correlations in growing and dividing cells. Within this framework, we derive an analytical expression for mRNA- protein correlations and show that explicit incorporation of growth and division substantially alters these correlations. The resulting relation is invariant to upstream transcriptional dynamics, and we validate it using stochastic simulations across multiple gene-regulatory architectures. Finally, we show that the derived predictions are consistent with the E. coli data.

Variability in gene expression among single cells and growing cell populations can arise from the stochastic nature of protein synthesis, which often occurs in random bursts. This study investigates the variability in the expression of a growth-sustaining protein, whose concentration is regulated by a negative feedback loop due to cell growth-induced dilution. We model the distribution of protein concentration using a Chapman-Kolmogorov equation for single cells and a population balance equation for growing cell populations. For single cells, we derive an explicit solution for the protein concentration distribution in state space and represent it as a Bessel function in Laplace space. For growing populations, we find that the distribution satisfies a Heun differential equation with singular boundary conditions. By addressing the central connection problem for the Heun equation, we quantify the population-level protein distribution and determine the Mathusian parameter, which characterizes population growth. This work provides a comprehensive analytical framework for understanding how stochastic protein synthesis impacts gene expression variability and population dynamics.

Successful development of multicellular organisms requires cells to transition between states with precise timing. Distinct cell states are often understood as being maintained by stabilizing regulatory networks, such that a complete cell-state transition requires network rewiring through partial dismantling of the current state and concurrent reconfiguration into a new one. Empirically, these transitions are often investigated by quantifying the gain or loss of expression of a small number of state-specific markers, frequently a single proxy. A general quantitative framework for comparing the kinetics of such transitions across experimental conditions is lacking. Here, we show that the delayed Weibull distribution provides a natural description of cell-state transition kinetics when transition is viewed as the cumulative consequence of many molecular events, whose timing may vary between cells and conditions, analogous to system failure in reliability theory. This formulation yields a compact model with interpretable parameters describing the delay before transition onset, the characteristic timescale of transition, the temporal form of the transition hazard, and the fraction of cells competent to respond. Together, this framework provides a practical and interpretable approach for quantifying the kinetics of cell-state transitions and how they are altered by perturbation.

A major challenge in neuroscience is to predict how currents in nanodomains affect voltage and ionic concentrations. Cable and Rall theory provide analytic current-voltage relations by neglecting concentration gradients, and the impact of concentration gradients is usually studied numerically with the Poisson-Nernst-Planck (PNP) model. A precise quantitative understanding of the combined dynamics remains limited because analytic current-voltage-concentration relations are missing. In this work we derive such relations using a novel approach based on cross-diffusion equations. For narrow cylindrical domains, we derive time-dependent and steady-state expressions that explicitly show how currents affect voltage and ionic concentrations. We find that the influx of only one ion can significantly change the concentrations of all the other ions even if no channels for these ions are present. After a current injection we compute a biphasic voltage transient where the small-time asymptotic corresponds to the steady-state solution of the cable equation. We show that the accuracy of cable theory prediction for the voltage depends on how the current is distributed among the various ions. Finally, we develop an iterative method to accurately compute steady-state profiles for voltage and concentrations using first-order results by subdividing a cylinder into small segments.

Phase separated condensates are recognized as a ubiquitous mechanism of spatial organization in cell biology. Biophysical modeling of condensates provides critical insights into the dynamics and functions of these subcellular structures that are difficult to extract via experiments. Here we present an efficient computational pipeline, CASPULE (Condensate Analysis of Sticker Spacer Polymers Using the LAMMPS Engine), to simulate and analyze the biological condensates made of sticker-spacer polymers. CASPULE implements a unique force field that combines traditional Langevin dynamics with a "detailed balance proof" protocol for single-valent bond formation between stickers. This framework allows us to study the non-trivial biophysics that emerge out of the single-valent sticker interactions coupled with the effect of separation in energetic contribution by stickers and spacers. We provide detailed documentation on how to setup the simulation environment, perform simulations and analyze the results. Through case studies, we highlight the utility and efficacy of our pipeline. Importantly, we provide statistical parameters to characterize the cluster size distribution often observed in biological systems. We envision this tool to be broadly useful in decoding the interplay of kinetics and thermodynamics underlying the formation and function of biological condensates.

The light absorption of the protein cryptochrome and its chromophore FAD is important for the regulation of circadian rhythms and in some species for sensing magnetic fields. To compute the absorption spectrum of chromophore, typically only a small region is treated quantum-mechanically due the high computational cost of spectroscopic calculations. We present a formalism that allows a quantum-mechanical treatment of not only the chromophore but also the neighbouring amino acids which differ from species to species. This is achieved by using the real-time time-dependent Hartree-Fock method. This method allows extending the quantum domain from typically only a few dozen atoms up to around 1,200 atoms for the largest calculations. The presented framework allows the treatment of neighbouring tryptophan residues or the cofactor molecule MTHF in the same calculation and allows to extract information of which regions absorb light depending on wavelength. The presented results also show that the environment around the chromophore FAD amplifies the light absorption in cryptochrome.

Alzheimers disease (AD) is characterized by the accumulation of amyloid-{beta} (A{beta}), yet the specific link between plaque burden and cognitive decline remains a subject of intense investigation. This paper presents a mathematical model that simulates the coupled dynamics of A{beta} monomers, soluble oligomers, and fibrillar species in the brain tissue. By modifying existing moment equations to include a dedicated conservation equation for A{beta} monomers, the model explores how various microscopic processes, such as primary nucleation, surface-catalyzed secondary nucleation, fibril elongation, and fragmentation, contribute to macroscopic disease progression. Central to this study is the concept of "accumulated neurotoxicity" as a surrogate marker of biological age, defined as the time-integrated concentration of soluble A{beta} oligomers. Unlike plaque burden, accumulated neurotoxicity cannot be reversed, and the harm it causes depends critically on the sequence of events that produced it. Numerical results demonstrate that while plaque burden and neurotoxicity both increase over time, their relationship is non-linear and highly sensitive to the efficiency of protein degradation machinery. Specifically, impaired degradation leads to a rapid advancement of biological age relative to calendar age. The model further identifies oligomer dissociation and fibril fragmentation as potential protective mechanisms that can counterintuitively reduce neurotoxic burden by diverting monomers away from the soluble oligomer pool. These findings provide a quantitative framework for understanding why individuals with similar plaque burdens may experience vastly different cognitive outcomes, underscoring the importance of targeting soluble oligomers early in therapeutic interventions.